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Abstract 



The calculation of the symmetry factor corresponding to a given Feynman 
diagram is well known to be a tedious problem. We have derived a simple 
formula for these symmetry factors. Our formula works for any diagram in 
scalar theory {(j)^ and 0^ interactions), spinor QED, scalar QED, or QCD. 
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I. INTRODUCTION 



In field theoretic calculations, amplitudes are calculated from the generating functional 
by performing functional derivatives. The Feynman rules are a handy short cut. One draws 
all the topologically distinct diagrams and uses a set of rules, derived from the generating 
functional, to represent the diagram as a integral over momentum variables. The Feynman 
rules tell us what factor to include for each propagator, vertex and loop momentum. How- 
ever, the overall numerical factor has to be determined separately. This factor is called the 
symmetry factor of the diagram One attempts to choose the coefficients of the interac- 
tion terms in the Lagrangian so that the symmetry factor is one. However, for any theory 
that has an interaction which contains more than one factor of the same field, there will be 
diagrams for which the symmetry factor is not one. Spinor QED is an example of a theory 
for which all symmetry factors are one. The interaction part of the Lagrangian has the 
form ijj'y^A^i/j. In this expression the electron, positron and photon fields appear once each. 
In contrast, in scalar QED, there is an interaction of the form (f)''A^A^(f). The presence of 
two powers of the photon field will result in non-unity symmetry factors for some diagrams. 
Similarly, 0^ theory, (p'^ theory and QCD will have non-unity symmetry factors for some 
diagrams. 



II. SYMMETRY FACTORS FOR SCALAR THEORIES 

We will write the symmetry factor as 

S = Af/V (1) 

where A/" and V stand for numerator and denominator respectively. We describe below how 
to determine the numerical values of Af and V. For simplicity, we begin with scalar theory. 
Consider the generating function for connected n-point functions for a scalar theory. We 
have, 



int 



i 6J{x) 



(2) 



W[J] = exp J J dy J{x)D{x - y)J{y) + i J dxC ^ 
We consider a theory with both cubic and quartic interactions: 

An. = + (3) 

The factors 3! and 4! in the expression above are included so that as many diagrams as 
possible will have a symmetry factor of one. We adopt the following short hand notation: 
Dxy := D{x — y); Jx := J{x); DJ^ ■= D^yJy] and JDJ := JxD^yJy The integrals over the 
spatial variables will be suppressed. In addition, we will drop all factors of i and not keep 
track of signs, since we are only interested in the numerical factor. 

We look at the calculation of the amplitude corresponding to the diagram in Fig. [1]. 
We obtain, 
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G{w,x,y,z) ^ ^^^^W[J]\j=o (4) 

OJwOJxOJyOJz 

where the notation ~ indicates that signs and factors of i have been dropped. We expand 
the exponential and keep only the terms that give the diagram we are considering. We 
obtain: 



5 5 5 5 f 5 V f 5 V f 5 




5J,,, 5JrS.L 5Jz [Sjaj \5Jh} [Sjrj ^ ^ 



The factor I) is given by: 



- = -(-)(- 

V 3! V4!y V2 



(6) 



This expression is obtained as follows: 

a) ^ from expanding exp(£j„t) = J2n ^(-Cmi)" and taking the third term in the expansion 
which contributes three factors of the interaction Lagrangian 

b) 3 factors of 1/4! from Cint (Eqn. (|^)) 

c) |y from expanding exp (^^ J dx J dy J{x)D{x — y)J{y)^ and keeping the term that 
contains eight propagators 

The factor Af is the number of ways that the functional derivatives can act on the factor 
[JDJY to produce the combination of propagators that corresponds to the diagram we 
are looking at. The first step is to include a factor 3! which comes from the possible 
permutations of the indices {a, h, c}. Including this factor allows us to restrict ourselves to 
the product of progagators: D axD awD ab{,Dhc)^ D acDhyD cz- Now we begin taking derivatives. 
We start with the derivatives that correspond to external legs. We obtain, 

16 ■ 14 ■ 12 ■ 10 ■ DJ^ ■ DJx ■ DJy ■ DJz ■ {JDjf (7) 

Note that there are no terms of the form D^w etc. since propagators of this form are not 
part of the diagram we are interested in. Proceeding with the first jj- we obtain, 

Da^ ■ DJx ■ DJy ■ DJz ■ {JDjf + Dax ■ DJ^ ■ DJy ■ DJz ■ {JDjf 

+8 ■ DJ^ ■ DJx ■ DJy ■ DJz ■ {JDjf ■ DJa (8) 

where we have differentiated -DJ^, -DJ^, and {JDJY in the first, second, and third expres- 
sions respectively. In the last expression, the 8 comes from acting on the {JDJ^. We 

■5 



proceed to the second t4-. We obtain 



Da^ ■ Dax ■ DJy ■ DJz ■ {JDjf 

+8 ■ Dau, ■ DJx ■ DJy ■ DJz ■ {JDjf ■ DJa 

+Dau, ■ Dax ■ DJy ■ DJz ■ {JDjf 

+8 ■ Dax ■ DJ^ ■ DJy ■ DJz ■ {JDjf ■ DJa 



+8 • Day, ■ DJx ■ DJy ■ DJz ■ {JDjf ■ DJ, 



+8 ■ Dax ■ DJy, ■ DJy ■ DJz ■ {JDjf ■ DJa 

+8 ■ 6 • DJy, ■ DJx ■ DJy ■ DJz ■ {JDjf ■ {DJaf (9) 
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and collecting like terms gives 

Da. ■ DJy ■ DJ, ■ {JDjf 
16 ■ ■ DJ, ■ DJy ■ DJ, ■ {JDjf ■ DJa 
16 ■ Da. ■ DJ^ ■ DJy ■ DJ, ■ {JDjf ■ DJa 
48 ■ DJ^ ■ DJ. ■ DJy ■ DJ, ■ {JDjf ■ {DJaf 

(10) 

There are two more ti-'s. The result is: 

{2 ■ S ■ Q + IQ ■ Q ■ 2 + IQ ■ Q ■ 2 + AS ■ 2) Dan, ■ Da. ■ DJy ■ DJ, ■ {JDjf ■ {DJaf (11) 
Combining, the coeffiecient obtained from the a derivatives is 

(2 ■ 8 ■ 6 + 16 ■ 6 ■ 2 + 16 ■ 6 ■ 2 + 48 ■ 2) = 12 • 8 ■ 6. (12) 



One must proceed in the same fashion with the 6, and c derivatives. The result of this 
lengthy procedure is: 

TV = 3! ■ (16 ■ 14 • 12 ■ 10) ■ (12 ■ 8 ■ 6) ■ (24 ■ 4 ■ 2) ■ 24 (13) 

Combining Eqns. (|I]), (H) and (|l^) we obtain, 

3! ■ (16 ■ 14 ■ 12 ■ 10) ■ (12 ■ 8 ■ 6) ■ (24 ■ 4 ■ 2) ■ 24 _ 1 

3! ■ (4!)3 . (2)8 ■ 8! ~2 ^ ^ 



The procedure described above, although straightforward, is incredibly tedious. The 
result can be obtained immediately using the formula: 

where: 

R = the number of ways to permute the internal indices and produce an identical set of 
propagators 

Di = the number pairs of progatators of the form D^a 

D2 = the numbers of pairs of propagators of the form -D^„ plus the number of factors of 
the form Dl^ (note: a factor D^^ contributes 1 to Di and 2 to - these propagators are 
counted twice) 

D3 = the number of triples of propagators of the form -Dj^„ 
D4 = the number of sets of propagators of the form D^^ 

We explain this notation by giving some examples. Figs. [2a-2c] show three diagrams 
with non-unity values of R. We list the R values and the indices or sets of indices that can be 
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permuted without producing a different set of propagators. Figs. [3a-3d] contain examples 
of diagrams witli different values of Di, D2, -D3 and D4. The values are listed beneath each 
diagram. Finally, we note that the factor (1/2)^^ in (^) is only different from one for the 
diagram in Fig. [4] . 

For the diagram in Fig. [1], the values of the parameters defined above are: R = 1, 
Di = = D4, = and D2 = 1 (because there is a pair of propagators of the form D^^). 



Substitution into ([T5|) gives 



in agreement with (|l^. Some more complicated diagrams are shown in Fig. [5]. 



III. GENERALIZATION TO OTHER THEORIES 

The formula ([T^) is valid for any diagram in (p^ theory, 0^ theory, spinor QED, scalar 
QED, or QCD. In this section we discuss the reason that the same formula works for different 
theories. It is not at all clear that this should be the case, since QED and QCD involve 
more than one type of field, including fields that are not self conjugate. Examples are the 
fermion fields in spinor QED {ip, ip) and the charged scalars in scalar QED (0, 0^). To 
handle fields of this type the generating functional (0) must be modified. We introduce 
separate sources for each field. A field that is not self conjugate is treated as two fields, and 
requires two sources (one for the field and one for its conjugate). For example, for QED 
the generating functional has the form (suppressing all indices, since they are not important 
when calculating symmetry factors): 



W[J] = exp\ I dx I dy [ -J^D^yJy + iK^G^yKy ] + i j dx Cint 



15 15 15 



i 5Jr i 5Kr i 5Kr 



(17) 



where K is the source for the field ip or (f)'^ , K is the source for the field ip or 0, and J is 



the source for the photon field. The inverse photon propagator is D{x 
fermion or scalar propagator is G{x — y). 



y) and the inverse 



A. Propagators 



First we look at propagators. Consider what happens when we take functional derivatives 
of the generating functional (p!7|). Compare the two expressions: 



5Jy 5Jx 

J 5_ 

5Ky 5K, 



K \ 



2 

■ G. 



xy 



D 



yx) 



D., 



xy 



xy 



(18) 



Since the photon is self conjugate, we have Dr^y = Dy^. There is only one source function for 
the photon field and the factor of 1/2 in front of the photon inverse propagator is necessary 
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to remove the factor of 2 that arises from taking the functional derivatives in either order. 
In contrast, the non-self-conjugate field has two independent degrees of freedom and thus 
we need two source functions, and no factor of 1/2 in front of the inverse propagator. 
The propagator has a direction that corresponds to the flow of the associated quantum 
number and thus does not satisfy Gxy — Gyx- The point is that, for both types of fields, 
the generating functional is constructed so that the propagator appears without spurious 
factors of two. 



B. Vertices 

In order to understand the behaviour of vertices in QED and QCD we look at a few 
examples. Consider the interaction between two photons and two charged scalars in scalar 
QED. The interaction Lagrangian has the form, 

JC^nt = -e^g^'^cp'^A^A^ (19) 

The corresponding vertex can be obtained by functional differentiation of the generating 
functional in the usual way. The result is V^^ = 2ie^g^„. The factor of two is a consequence 
of the fact that the two photon fields are inter changable. Now consider the calculation of 
the symmetry factor. Niavely we would write the vertex in functional derivative form as 

Consider what happens if we proceed with the calcuation of the symmetry factor following 
the example for scalar fields described in section II. For a scalar theory, the factor of two 
which appears in the vertex (from the interchange of two fields), would have been included 
in the factor J\f. Thus, if we want to treat the photons as scalars (for the purposes of 
calculating the symmetry factor) we have to include an extra factor 1/2 in the interaction. 
The 'scalarized' interaction has the form: 




^ _d__S_ 

JkJk 



(21) 



Now, consider the three gluon vertex in QCD. This vertex is obtained from the part of 
the interaction Lagrangian that contains three gluon fields which can be written: 

Cint = -gfabcd^AlA'^A''^ (22) 

The corresponding vertex is: 

rJp(P> ^) = 9fabc[g,,v{p - q)p + gup{q - r)^ + gpt,{r - p)^] (23) 

Note that the vertex contains six terms. In a scalar theory, these six terms would be identical 
and would give rise to a factor of six, were it not for the fact that the interaction Lagrangian 
contains an explicit factor of 1/3! which exactly cancels this factor of six. Therefore, to 
obtain the symmetry factor, we can use a 'scalarized' interaction Lagrangian of the form 
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where J is the source for the gluon field, and all indices are ignored. The point is that this 
interaction has exactly the same form as the interaction for 0^ theory considered above. 
The four gluon vertex works the same way: The interaction Lagrangian is 

g'i y \ r 

C-int = -^fabcfadeQiiugXT^b^^c^d^^e (25) 

The Feynman rule for the vertex contains six terms. Combining the factor of 1/4 in the 
interaction Lagrangian and a factor of 1/6 so that we can treat the gluons as scalars, we 
need a factor of 24 = 4! in the 'scalarized' form of the interaction vertex: 

This interaction has the same form as the interaction for 0'^ theory considered above. 

All interactions in spinor and scalar QED and QCD can be written in a 'scalarized' form 
by including a factor 

^i—, (27) 
where Ni is the number of fields of type i in the vertex. 



C. Symmetry factors 

Using the discussion in the sections above about the behaviour of propagators and ver- 
tices for non-scalar theories, we can consider taking functional derivatives to construct a 
given diagram, following the method outlined in section II. It is straightforward to see that 



symmetry factors are given by the same formula as for scalar theory (15) by using the rules 
below (in all figures, solid lines are uncharged scalars, solid lines with arrows are fermions 
or charged scalars (or ghosts), and dotted lines are photons or gluons): 

a) Factors of G make no contribution to the factors Di, D4 in ([T5|). For example, 
the first diagram in Fig. [6] has D2 = 1 and thus S = 1/2, while the second diagram in Fig. 
[6] has ^2 = and 5 = 1. 

b) The fact that Gxy 7^ Gyx must be taken into account when calculating R factors. For 
example, the first diagram in Fig. [7] has R = 2, but the second diagram in Fig. [7] has 
R=l. 

Fig. [8] gives some more examples of diagrams for which the symmetry factor can be 
calculated using ([T5|). 
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IV. CONCLUSION 

The calculation of symmetry factors is a straightforward but tedious business. The 
formula derived in this paper ([T5| ) provides a simple way to obtain the symmetry factor for 
any diagram in 0^ theory, 0^ theory, spinor QED, scalar QED, or QCD. 
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Fig. [4]: R=l, Di = l, D2 = 2, 03 = 0,0^ = 




(a) (b) (C) (d) 




Fig. [5]: Some diagrams for scalar theories. 
Symmetry factors are: a) 1/6; b) 1/12; c) 1/2; d) 1/2; c) 1/6; f) 1/12; g) 1/4; h) 1/192 
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Fig. [6] 



D factors for theories with non-self-conjugate fields 




Fig. [7]: R factors for theories with non-self-conjugate fields 




Fig. [8]: Some diagrams with non-self-conjugate fields. 
Symmetry factors are: a) 1; b) 1/2; c) 1/2; d) 1; e) 1/4; f) 1 
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